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INTRODUCTION 
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The  interaction  between  the  chemical  species  which  make  up  the 
upper  atmosphere  and  the  energy  sources  which  drive  the  atmospheric 
chemistry  are  primarily  two-body  2nd  order  reactions  between  the  resident 
chemical  species  and  the  flux  of  incoming  photons  and  energetic  particles. 
The  rate  constants  for  the  reactions  are  expressed  as  interaction  cross- 
sections  which  can  and  often  do  vary  rapidly  with  available  energy  in 
addition  to  being  different  for  each  chemical  species.  Treating  these 
reactions  in  a straightforward  way  as  two-body  reactions  and  adding  them 
to  any  atmospheric  reaction  set  would  result  in  an  essentially  uncalcu- 
latable  set.  The  deliberately  broad  band  source  models  described  in 
this  report,  which  model  energy  inputs  with  only  47  solar  bands  plus 
photoelectrons,  precipitated  electrons,  and  galactic  cosmic  rays,  would 
be  prohibitively  unwieldy  if  combined  with  the  45  ^asic  chemical  species- 
photon  interaction  reactions  in  th^  ECLIPSE-66  set,  or  the  60  basic 
reactions  in  the  BENCHMARK-76  set.  After  allowing  for  those  rate 
constants  which  are  zero  in  the  present  sets,  or  would  be  zero  in  the 
full  set,  we  get  an  increase  from  perhaps  40  basic  reactions  to  about 
2000  for  these  two  cases. 

It  is  therefore  customary  to  reduce  these  photon-species  reactions 
to  psuedo-first  order  reactions  by  taking  the  sum 

k = y <)>.o. 
n ‘ Ti  in 

l 

where 

n = the  number  of  the  pseudo  first  order  reaction 

i 3 the  wavelength  band  number 

k 3 psuedo-first  order  rate  constant 

4>  3 photon  flux 

a 3 interaction  cross  section. 


*F.  E.  Niles  and  I.  L.  Chidsey,  "Comparison  of  Measured  and  Calcu- 
lated D-Region  Electron  Density  for  Solar  Eclipse  of  12  November  1966," 
and  I.  L.  Chidsey  and  F.  E.  Niles,  "Comparison  of  Measured  and  Cal- 
culated D-Region  Positive-Ion  Number  Density  for  Solar  Eclipse  of 
12  November  1966,"  Trans.  AGU  56,  997  (1975). 

2 

J.  M.  Heimerl  and  F.  E.  Niles,  "BENCHMARK-76  Model  Computations  for 
Disturbed  Atmospheric  Conditions,  I.  Input  Parameters,"  BRL  Report 
No.  2022,  October  1977;  "BENCHMARK-76  Model  Computations  for  Disturbed 
Atmospheric  Conditions,  II.  Results  for  the  Stratosphere  and  Mesosphere," 
BRL  Report  ARBRL-TR-02050,  March  1978;  and  "BENCHMARK-76  Model 
Computations  for  Disturbed  Atmospheric  Conditions,  III.  Results  for 
Selected  Excitation  Parameters  at  60  km,"  BRL  Report  ARBRL-TR-02051 , 

April  1978.  (AD  Nos.  A050355,  A054325,  and  A054376) 
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This  report  documents  how  we  obtained  these  psucdo-first  order  rate 
constants  for  two  studies,  ECLIPSE-66  and  BENCHMARK-76. 

For  the  ECLIPSE-66  study  the  time  variation  of  these  reaction  rate 
constants  was  critical,  and  they  were  computed  by  a subroutine  package 
which  was  a part  of  the  ECLIPSE-66  computer  program. 

For  the  BENCHMARK-76  study  the  time  variations  were  not  of  interest 
and  the  subroutine  package  was  used  as  a separate  program.  The  output 
of  this  "stand-alone"  program  was  a table  of  psuedo-first  order  rate 
constants  as  a function  of  altitude  and  is  recorded  as  Table  9 of  the 
Appendix.  [See  also  Table  3 of  reference  2 (I).] 


ENERGY  SOURCES 

The  code  contains  models  of  several  sources  of  energy:  the  sun, 
including  x-rays,  and  ultraviolet  and  visible  light;  precipitated  electrons; 
and  galactic  cosmic  rays.  The  solar  model  requires  a separate  atmospheric 
attenuation  model;  the  precipitating  electron  model  includes  atmospheric 
attenuation;  the  cosmic  ray  model  ignores  it.  The  eclipse  study  also 
required  an  eclipse  model. 


NON- ION I ZING  SOURCES 

The  only  source  of  non-ionizing  energy  in  the  code  is  the  sun. 
Although  the  visible  light  and  most  of  the  ultraviolet  light  are  unable 
to  produce  the  ionization  we  are  trying  to  study,  the  enormous  amounts 
of  non-ionizing  radiation  relative  to  the  ionizing  radiation  drive  most 
of  the  neutral  chemistry  and  can  strongly  modulate  the  charged  particle 
chemistry. 

The  Solar  Model,  which  is  given  in  Table  1 of  the  Appendix  is  based 
on  that  used  in  the  OPTIR  II  code.  The  OPTIR  II  Solar  ModeJ  covers 
the  ultraviolet  from  100  to  |00  nm  in  30  bands  each  0.2S  ym”  wide. 

(The  inverse  micrometer  (ym  ) is  a measure  of  wave  number  or  energy 
similar  to  the  inverse  centimeter.)  We  chose  to  base  our  solar  model 
on  the  OPTIR  II  model  because  it  is  made  up  of  a few  wide  bands  rather 
than  many  narrow  bands  and  is  given  in  units  of  photon  flux  rather  than 
of  energy  flux.  This  model  was  extended  on  the  short  wavelength  side 
by  ten  x-ray  bands:  three  bands  0.2  nm  wide  from  0.4  to  1.0  nm,  and 
seven  bands  nominally  1.0  nm  wide  from  3.1  to  10.0  nm.  The  x-ray  photon 
fluxes  for  bands  1-3  were  taken  from  Figure  3 of  Sears  review  of  November 


^0.  P.  Manley,  H.  J.  P.  Smith,  Y.  M.  Treve,  J.  W.  Carpenter,  T.  C. 
Degges,  and  L.  R.  Doone,  "OPTIR  II,"  AFCRL-71-0528  (1),  AFCRL, 
Bedford,  MA,  1971. 
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I960  eclipsq  measurements  with  the  conversion  from  ergs  cm  “ sec  to 
photons  cm  “ sec  computed  for  the  center  of  each  band.  Because  Scars' 
energy  flaxes  for  bands  4-10  are  approximately  5 times  those  recommended 
by  Swider  for  'quiet'  solar  conditions  we  used  5 times?Swidey's  'quiet' 
modql  for^our  model.  Again,  we  converted  from  ergs  cm  “ sec  to  photons 
cm  “ sec  at  the  band  center. 

The  long  wavelength  limit  of  the  OPTIR  solar  model  was  extended 
from  400  nm  to  735  nm  by  adding  seven  bands,  the  width  of  each  being 
about  10%  of  cegter  wavelength.  These  were  obtained  by  combining  bands 
from  Ackerman's  model,  as  shown  in  Table  1. 

Our  47  band  model  requires  about  2200  words  of  memory  when  used 
with  our  5 absorbing  species  and  60  photolytic  reactions.  Direct  use 
of  Ackerman's  model  would  have  quadrupled  both  memory  requirements  and 
the  number  of  photon  interaction  calculations  such  as  absorption,  dis- 
sociation, etc. 


IONIZING  RADIATIONS 

For  computational  purposes  it  was  convenient  to  consider  four  sources 
of  ionization:  solar  ultraviolet,  solar  x-rays,  galactic  cosmic  rays, 
and  precipitating  electrons.  The  solar  ultraviolet  ionization  source 
comprises  bands  11-20  of  the  solar  model  mentioned  above.  The  solar 
x-ray  source  comprises  bands  1-10  of  the  solar  model  and  includes  ioniza- 
tion by  two  processes:  primary  photo-ionization  and  secondary  ioniza- 
tion due  to  energetic  photo-electrons. 

The  Photoelectron  Pair  Production  Model  supplements  the  x-ray 
portion  of  the  solar  model.  The  x-ray  photons  have  considerable  energy 
beyond  that  required  to  ionize  atmospheric  chemical  species,  and  this 
energy  appears  as  kinetic  energy  of  t{ie  photoelectrons  plus  a small  con- 
tribution from  multiply-charged  ions.  These  photoelectrons  plus  any 
multiply-charged  ions  create  secondary  ionization  in  the  vicinity  of  the 
primary  photon  absorption.  This  secondary  ion-pair  production  is 
given  by 


4R.  D.  Sears,  "Analysis  of  the  1966  Solar  Eclipse  Data,"  LMSC  D246526, 
Lockheed  Palo  Alto  Research  Lab,  Palo  Alto,  CA,  1972. 

bW.  Swider,  Jr.,  "Ionization  Rates  Due  to  the  Attenuation  of  1-100A 
Nonflare  Solar  X-Rays,"  Rev.  Geophys.  7_,  No.  3,  573-594  (August  1969). 

Si.  Ackerman,  "Ultraviolet  Solar  Radiation  Related  to  Mesospheric 
Processes,"  in  Mesospheric  Models  and  Related  Experiments,  G.  Fiacco, 
ed..  Springer  Verlag,  N.Y.,  149-159  (1971). 
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TABU-  1.  DERIVATION  OF  URL  BANDS  41  - 47  FROM  ACKERMAN'S  SOLAR  MODEL6 


Ackerman's 

Bands 

BRL 

Band 

Photons 

,n16,  -2  -K 

10  (cm  sec  ) 

Wavelength 

Range  (nm) 

113-118 

41 

1.145 

400-435 

119-126 

42 

1.851 

435-470 

127-134 

43 

2.015 

470-510 

135-144 

44 

2.581 

510-560 

145-156 

45 

3.237 

560-620 

157-165 

46 

3.040 

620-680 

166-171 

47 

3.528 

680-735 

where 


Pairs0l  = 0.03  (Y  E . a.  I.)  - 1 
Phe  1 i i J 


PairSp,  ^ = photo-electron  pair  production 


E. 

l 


o . 

l 


I. 

l 

0.03 


Phe 

= energy  per  photon  in  electron  volts 

-22  2 

= ionization  cross  section  in  10  M 

-2  -1 

= the  attenuated  photon  flux  in  photons  cm  sec 

e the  approximate  ionization  efficiency  in  ion-pairs 
per  electron  volt. 

i = the  solar  x-ray  band  number. 

The  one  in  this  expression  accounts  for  the  initial  ionization  reaction. 

- 1 8 

The  Galactic  Cosmic  Ray  ionization  constant  used  was  3.6  x 10 
-1  -3  -1 

molecules  cm  sec  and  was  obtained  by  linear  interpolation  between 
Sears'  equations  1 and  2 

Q 


<GCR  max 


= 4(30  + 250  sin4*)  x 10'2°  n 


Q„„n  . = 4(30  + 1000  sin4*)  x 10"20  n 

XGCR  min  v J 

where  * = the  geomagnetic  latitude 

n = number  density  per  cm^ 

for  a geomagnetic  latitude  of  32.4°S4  and  is  consistent  with  Swi^er's^ 
table  7d.  A parameterization  of  cojgic  ray  ionijati^n  variation  gives 
an  ionization  ^onstant  of  3.0  x 10  molecules'  cmJ  sec'  . The  latitude 
given  in  Sears  as  the  geomagnetic  latitude  appropriate  for  the  eclipse 
data  used  is  apparently  the  geographic  latitude,  and  the  appropriate 
geomagnetic  latitude  is  17.7*  South.  Thus,  the  value  used  for  the 
Galactic  Cosmic  Ray  ionization  constant  is  probably  high  for  the  ECLIPSE 
and  low  for  the  BENCHMARK  runs.  However,  the  data  tabulated  by  Swider 
indicate  that  this  value  of  the  ionization  constant  is  not  unreasonable. 


The  cosmic  ray  model  makes  no  allowance  for  attenuation  of  the 
cosmic  ray  flux  by  the  atmosphere,  and  therefore  becomes  increasingly 
inaccurate  below  30  km. 

^M.  G.  Heaps,  "The  Ion  Pair  Production  Due  to  Cosmic  Rays,"  BRL 
Memorandum  Report  2641,  June  1976.  (AD  #A027629) 

g 

D.  C.  Baker,  "Ionospheric  D-Region  Parameters  from  Blunt  Probe 
Measurements  during  a Solar  Eclipse,"  Scientific  Report  #334,  pp.  45-46, 
Penn.  State  University,  February  1,  1969. 
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The  Precipitating  Electron  pair  production  model  is  based  on  that 
given  by  Sears.4  In  his  appendix  B,  which  is  reproduced  in  Table  3 of 
the  appendix,  he  gives  the  average  pair  production  rate  versus  altitude 
for  the  South  Atlantic  anomaly  at  eclipse  time.  In  his  Figure  1 he  also 
plots  the  maximum  pair  production  rate,  which  is  about  12  times  greater. 
Both  of  these  rates  are  based  on  extrapolations  from  trapped  electron 
flux  measurements  made  on  September  30,  1966,  some  six  weeks  before  the 
eclipse. 

Scars'  table  gives  data  from  54  to  100  km  only  and  the  model  was 
extended  above  100  kin  by  assuming  that  the  pair  production  fell  with  a 
constant  scale  height  above  98  km.  The  model  was  extended  below  54  km 
by  assuming  that  the  scale  height  decreased  at  a constant  rate  below 
58  km. 


ECLIPSE  FUNCTIONS 

One  additional  complication  arises  for  the  eclipse.  The  variation 
in  energy  input  must  be  known  if  we  are  to  accurately  model  the  effects. 
Sears  gives  several  eclipse  functions:  for  visible  light  he  gives  a 
table  of  values  computed  from  the  expression  for  a uniformly  bright 
circular  disc  eclipsed  by  a disc  of  the  same  subtended  angle  - 

Q = Qo  2 (arcsin  (t/d)  + (t/d)  (1  - (t/d) 2) 1/2)/n 


where 


Qo  = the  uneclipsed  flux 
t = time  before  or  after  totality 
d = the  semiduration  of  the  eclipse. 

9 

This  eclipse  function  is  strictly  true  only  in  the  far  infrared. 

For  most  wavelengths  an  optical  depth  of  unity  occurs  below  the  tempera- 
ture minimum  in  the  solar  atmosphere.  Near  the  center  of  the  solar  disc 
the  ray  patljig  are  vertical  and  the  spectrum  resembles  that  of  a 6700  K 
black  body.  Near  the  limb  the  ray  paths  are  oblique,  unity  optical 
depth  occurs  at  a higher,  cooler  ^titude,  and  the  effective  radiation 
temperature  is  5400K  at  the  limb.  Sears'  ultraviolet  eclipse  function, 
which  reflects  this  limb  darkening  is  given  in  Table  6 of  the  Appendix. 
For  Lyman  alpha  the  effects  of  limb  darkening  are  counteracted  by  reso- 
nant scattering  from  the  hydrogen  in  the  solar  corona  and  the  geocorona 
and  Sears  recommends: 

9G.  Abetti,  The  Sun,  Macmillan,  N.Y.,  1951,  Fig.  85,  p.  282. 

10G.  Abetti,  ibid. , Fig.  86,  p.  285. 
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for  t < 3400  seconds 


Q,  = Q . + .001 

Ly  xv  i s 


= Q . + .002 

xv  1 s 


for  t < 1700  seconds. 


The  x-ray  flux  originates  to  a large  extent  in  the  solar  corona  and 
is  highly  variable  in  space  and  time.  Sears’  recommended  eclipse  function 
is  given  in  the  Appendix  (Table  6)  and  is  quite  asymmetric. 


During  totality  the  eclipse  functions  used  were 


Q 


visible 


Q 


ultraviolet 


0. 


Q 


Lyman  alpha 


(.001  + .001*  t/d')Q 


where  t = time  from  the  center  of  totality  and  d'  5 the  semiduration 
of  totality.  For  x-rays  we  interpolated  linearly  between  the  values  at 
2nd  contact  (the  start  of  totality)  and  3rd  contact  (the  end  of 
totality) . 

9 

We  were  strongly  influenced  by  Abetti’s  limb  darkening  curves 
when  assigning  our  solar  bands  to  these  eclipse  functions.  We  used  the 
ultraviolet  eclipse  function  for  the  near  ultraviolet  and  throughout 
the  visible  (our  bands  28-47,  174-740  nm) . We  used  the  Lyman  alpha 
function  for  our  bands  11-24,  100-148nm,  and  the  "visible"  eclipse 
function  only  for  bands  25-27,  148-174  nm. 


ATMOSPHERIC  ATTENUATION  COMPUTATIONS 

The  solar  model  gives  us  the  photon  flux  at  the  top  of  the  atmosphere, 
but  the  information  we  need  is  the  photon  flux  at  some  lower  level  where 
we  are  simulating  the  atmospheric  chemistry.  We  model  the  atmospheric 
attenuation  of  the  solar  flux  by  Beer's  Law 

T T -No 

I = I e 
o 

-2-1 

where  I = initial  photon  flux  in  photons  - cm  - sec 

N = line  of  sight  number  density  per  square  centimeter 

0 = absorption  cross  section  in  square  centimeters 

1 = the  attenuated  photon  flux  one  band  at  a time. 

Because  there  are  more  than  one  absorbing  species  we  replace  No 

by  j Ni  o i to  get 
i 
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Seven  atmospheric  species  enter  into  the  attenuation  calculations: 

Argon,  carbon  dioxide,  atomic  and  molecular  nitrogen,  atomic  and 
molecular  oxygen,  and  ozone.  The  absorption  cross  sections  given  in 
Table  1 of  the  Appendix  were  used  both  for  the  ECLIPSE  calculations  and 
for  t he j. BENCHMARK”  calculations.  The  x-ray  cross  sections  are  from 
Swider.'  The  cross  sections  for  the  other  wavelengths  can  be  found  by 
looking  for  the  same  species  in  the  lists  of  photolytic  reactions  in 
the  Appendix. 

The  species  number  densities  are  the  same  as  those  used  as  initial 
data  for  the  AIRCHEM1  and  BENCHMARK  computations,  “ but  the  absorbing 
species  density  model  has  been  expended  from  0-500  km  using  the  DNA 
reaction  rate  handbook  DNA  1948H.  In  some  cases  densities  have 
been  extrapolated  beyond  those  given  in  the  sources  to  avoid  sudden 
jumps  in  the  computed  flux  level  for  small  changes  in  ray  path. 

The  column  densities  of  the  absorbing  species  were  found  by  integra- 
ting the  column  densities  along  the  ray  path  from  its  lowest  point  until 
the  OPTIR  convergence  criterion  of  a less  than  1 percent  change  in  the 
column  density  per  integration  step  had  been  met  for  all  absorbing  species. 
The  scale  path  length  (i.e.,  the  projection  of  the  scale  height  on  the 
ray  path)  was  then  computed  for  each  absorbing  species  and  was  used  to 
extrapolate  the  column  density  to  the  top  of  the  atmosphere.  The  column 
densities  were  then  checked  again  to  ensure  that  either  the  extrapolated 
column  density  was  less  than  .01  gm  cm  (i.e.,  negligible)  or  that  the 
extrapolation  caused  an  increase  of  less  than  10  percent  in  the  column 
density.  If  any  species  failed  both  of  these  criteria  the  OPTIR  criterion 
was  halved  and  integration  along  the  ray  path  was  resumed. 

Two  problems  with  the  absorption  computations  must  now  be  considered. 


I . Time  Variation  of  Absorbing  Species 

In  general  the  chemistry  may  cause  the  number  density  of  an  absorb- 
ing species  to  change  with  time.  For  the  present  case,  only  one  absorb- 
ing species,  ozone,  changes  enough  to  cause  a problem.  If  the  ray  path 


^E.  L.  Lortie,  M.  D.  Kregel  and  F.  E.  Niles,  "AIRCHEM:  A Computational 
Technique  for  Modeling  the  Chemistry  of  the  Atmosphere,"  BRL  Report 
No.  1913,  August  1976  (AD&A030157) . 

12F.  E.  Niles  and  J.  M.  Heimerl,  "Selected  Neutral  Species  Profiles 
0-100  km,"  BRL  Memo  Report  No.  2767,  July  1977  (AD#A042620) . 

^DNA  Reaction  Rate  Handbook,  DNA  1948H,  Chapter  2,  Revision  1,  DASIAC, 
GE-TEMPO,  Santa  Barbara,  CA,  1972. 


14 


is  nearly  vortical  so  that  most  of  the  absorption  occurs  near  the 
altitude  of  the  chemical  computations,  and  the  species  in  question 
causes  only  a small  amount  of  absorption,  we  may  make  the  assumption 
that  the  variation  in  line  of  sight  density  is  proportional  to  the  j 
variation  in  computed  density.  This  .assumption  holds  for  the  hCLII’SI. 
the  ray  path  was  within  -10°  of  vertical  and  the  computations  were  fqr 
altitudes  > i>()  km,  well  above  the  ozone  maximum.  for  the  BENCHMARK' 
runs  the  density  of  the  absorbing  ozone  was  held  constant;  the  solar 
inputs  were  computed  once  for  an  overhead  sun  using  the  initial  model 
densities.  This  assumption  would  not  hold  for  ray  paths  which  dipped 
below  the  ozone  maxi  mum,  and  would  he  grossly  violated  for  sunrise/ 
sunset  studies. 

I 1 . Choice  of  Effect  i ve  tiros s Sect  ions 

The  relatively  wide  bands  in  our  solar  model  lead  to  another  problem. 
The  cross  sections  of  absorbing  species  can  vary  by  several  orders  of 
magnitude  over  one  of  our  bands,  making  our  choice  of  the  proper  effec- 
tive cross  section  difficult.  A particularly  severe  problem  is  presented 
by  the  Schuman-Runge  bands  of  molecular  oxygen,  which  roughly  corres- 
pond to  our  bands  27-31  (167-210  nm) . Oxygen  also  has  a strongly  varying 
cross  section  over  our  bands  13-17  and  19  (105-121  nm  and  125-129  nm). 

For  our  hand  18,  which  is  dominated  by  Lyman  alpha,  we  use  cross  sections 
for  Lyman  alpha.  Parts  of  this  discussion  also  apply  to  N,  absorption, 
bands  11-23  (100-148  nm) . 

Fortunately,  it  can  he  shown*'’’'*1  that  the  absorption  due  to  an 
isolated  gaussian  line  which  is  narrow  with  respect  to  the  measuring 
instrument  has  an  effective  absorption  cross  section  which  varies 
inversely  witn  the  square  root  of  the  absorber  density,  and  that  this 
relationship  holds  provided  the  opt ical ^ensity  at  the  peak  of  the  line 
is  at  least  unity.  It  can  also  be  shown  that  this  relationship  holds 
approximately  for  triangular  lines,  for  randomly  spaced  arrays  of  lines, 
and  for  evenly  spaced  arrays  of  lines.  £his  relationship  has  been 
experimentally  verified  by  Blake  at  al.  for  a range  of  absorption  suffi- 
cient to  vary  the  light  intensity  in  the  Schuman-Runge  bands  by  two  orders 
of  magnitude. 
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(1.  Kockarts,  "Penetration  of  Solar  Radiation  in  the  Schuman-Runge 

Bands  of  Molecular  Oxygen,"  in  Mesospheric  Models  and  Related  Experiment s , 

G.  Fiacco,  ed. , Springer  Verlag,  N.Y.,  pp.  160-176. 
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A.  J.  Blake,  J.  II.  Carver,  and  G.  N.  Haddad,  "Photoabsorption  Cross 
Sections  of  Molecular  Oxygen  Between  1250A  and  2350A,"  J.  Quant. 

Spect.  and  Rad.  Transfer  6,  453-454  (1966). 

*6J.  R.  Nielson,  V.  Thornton,  and  E.  B.  Dale,  "Absorption  Laws  for 
Gases  in  the  Infrared,"  Rev.  Mod.  Phys.  1_6-  307  (1944). 

* A.  E.  S.  Green  and  P.  .J.  Wyatt,  Atomic  and  Space  Physics, 

8-3  and  8-4,  pp.  410-429,  Addison-Wesley , N.Y.,  1968. 


Sect  ions 


Wo  now  consider  whnt  cross  sections  to  use  for  those  bands  for  which 

the  cross  section  varies  too  widely  for  one  cross  section  to  be  suitable. 

Wo  assume  that  the  unattenuated  photon  flux  is  relatively  constant  across 

the  band.  We  also  assume  that  there  is  some  reasonable  distribution  of 

cross  sections  between  a , the  maximum  cross  section  within  the  band, 

max 

and  the  minimum  cross  section.  Consider  dividing  a band  into  sub- 

bands  whose  boundaries  are  defined  by  o , o ...  a ..., 

^°min’  °min'  Wc  accept  a cross  section  distribution  as  reasonable 

if  the  sub-band  intervals  are  approximately  equal.  If  this  is  true 


n = 


— -r-  £n(o  /o  . ) 
1.4  max  min 


and  1.4  = Jin (4)  . 


Then  when  the  number  density  of  an  absorber  along  the  line  of  sight  is 
less  than  l/amax  only  the  peaks  are  contributing  to  the  absorption  and 


o 

e 


0.5  a 


max 


/n 


where  oe  = the  effective  cross  section,  i.e.,  the  one  we  should  use  in 

our  absorption  calculations.  In  other  words,  1/2  of  a is  effective 
r ’ max 

over  1/nth  of  our  band  or  (l/2n)th  of  a is  effective  over  the  whole 

, , max 

band. 


1 8 

Hudson  indicates  that  in  much  of  the  ultraviolet  region  the 
measured  cross  sections  are  instrumental  bandwidth  dependent,  and  that 
in  general  the  peak  cross  sections  should  he  greater  and  the  minimum 
cross  sections  should  be  lower.  We  will  therefore  assume  that  for  the 
most  intense  sub-band 


a=a  /n  = 1.4  o /J,n(o  / o . ). 

e max  max  v max  min 

If  the  number  density  is  such  that  l/o  <<  N <<  l/o  . and  if  the 

max  min 

t U 

inverse  square  law  does  hold,  then  oq  = oQ  NQ 2/N 2 where  NQ  is  the  line  of 

sight  density  when  o = o...  If  we  let  o.  = o when  N = l/o 

e 0 0 e max 


o =o 
e max 


1.4/ (£n(o  /o  . )(o  N)  2) 
y v max'  mm'  ^ max  ’ ’ 


- 1.4  o 2/(N2  tn(o  /o  . )) 


max  min' 
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R.  D.  Hudson,  " Review  of  UV  Cross  Sections 
Phys.  9,  No.  2,  304-406  (1971). 
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If  the  number  density  is  very  large  then  N >>  l/o  . and  o > o . . 

1 ’ b min  e min 

We  now  have  3 expressions;  it  would  be  simpler  if  we  only  had  one. 


Clearly,  we  can  combine  the  latter  two  cases  by  simple  addition  to 


get 


o = 1.4  (o  /N)  2/Hn(o  /o  . ) + a . 
e max  ' max  mm  min 


for  all  N >>  l/o 


max 


Similarly,  we  can  combine  the  first  two  cases  by  taking  the  har- 
monic sum  1/S  = 1/a  + 1/b. 


If  we  let 


a = 1 . 4 o /X,n(o  /o.) 

max  max  mm 

b = 1.4  o 2N  2/in(a  /o  . ) 
max  max  min 


and  S = o 


then 


o = 1.4  o /(£n(o  /o  . ) (o  " N 2 + 1)) 
e max  max  min'  max 


for  N < l/o  . . 

mm 

Since  o = a . for  N > l/o  . we  can  add  this  to  get 
e mm  mm 


j,  i , 

"2  w'2 


o=o.  ♦ 1.4  o /(Hn(o  /o  . ) (o  N + 1)) 
e mm  max  max  min  max 


for  all  absorption  calculations. 


Because  we  now  have  two  cross  sections,  o and  o . , rather  than 

max  mm 

one  for  some  bands  we  now  have  a minor  storage  and  retrieval  problem. 

For  O0,  and  only  for  those  bands  where  the  problem  existed,  o was 
z max 

stored  as  an  absorption  cross  section  and  o . as  a dissociation  cross 
r mm 

section.  For  N_,  a was  stored  as  an  absorption  cross  section  and 
* max_24  2 
a . was  set  to  3 x 10  cm  . 
mm 

For  CO2  there  is  a different  problem.  The  CO.,  absorption  cross 
section  decreases  rapidly  with  increasing  wavelength  across  the  0~ 
Schuman-Runge  bands.  The  0 2 absorption  is  greater  and  strongly  modu- 
lates the  solar  source  seen  by  the  CO^.  The  effective  C 0^  absorption 
is  therefore  also  modulated  by  the  absorption.  To  simulate  this 
modulation,  the  C0-,  cross  section  was  made  to  vary  by  the  same  pro- 
portion of  its  range  between  a and  a . that  the  CL  cross  section 
r b max  mm  2 

did.  Because  the  C0o  cross  section  was  decreasing  monotomically  with 
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for  band  n+1 


increasing  wavelength,  o . for  band  n could  be  set  to  o 

nun  max 

and  so  o was  stored  for  each  band, 
max 


Note  that  these  attempts  to  allow  for  square  root  low  absorption 
result  in  some  cross  sections  for  09,  N£,  and  C09  which  seem  to  hrve 
unusually  high  or  low  values. 

Because  of  an  algebraic  error  the  factor  1.4  in  the  final  expres- 
sion for  the  cross  section  was  replaced  by  a 2 in  the  computer  program. 
The  effect  of  this  error  will  be  discussed  in  the  section  on  photolytic 
reactions . 


Atmospheric  attenuation  of  the  precipitating  electron  flux  is 
included  in  the  model.  Atmospheric  attenuation  of  the  galactic  cosmic 
ray  flux  is  small  above  30  to  35  kilometers  and  the  effects  of  attenua- 
tion are  therefore  ignored. 


PHOTOLYTIC  REACTIONS 


These  are  reactions  which  are  driven  by  photons  and  which  decompose 
a chemical  substance.  The  rate  constants  of  the  photolytic  reactions 
are  computed  using  the  expression 


R = y a . 
n v nj 
J J 


I . 
J 


where 


n = reaction  number 


j = wavelength  band  number 


R = reaction  rate  constant 
n 


a . 
nj 


reaction  and  band  specific  cross  section 


I. 

J 


attenuated  photon  flux. 


The  reactions,  cross  sections  and  references  for  the  ECLIPSE  study  are 
listecj  in  Tables  4 and  5 of  the  Appendix,  and  those  for  the  BENCHMARK 
study  are  listed  in  Tables  7 and  8 of  the  Appendix. 


Because  the  various  ionizing  mechanisms  differ  mainly  in  their  energy 
sources  we  have  been  able  to  reduce  the  number  of  chemical  reactions. 
Ionization  by  Energetic  Particles  can  be  described  by  chemical  equations 
which  are  formally  the  same  as  the  equivalent  photo- ionization  reaction 
equations.  These  equivalent  reactions  were  combined  and  the  rate  con- 
stant for  each  combined  reaction  was  obtained  by  adding  those  of  its 
constituent  reactions.  For  example,  there  is  only  one  ionization 
reaction  for  molecular  oxygen 
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0,  + hv  ► 0,+  + c 

with 

R = U + R + R , + R + R 

n uv  x phe  prc  gcr 

where 

hv  photons  from  all  ionization  sources 
R reaction  rate  constant 
n - reaction  number 
uv  = ultraviolet 

x = x-rays  (primary  ionization  only) 
phe  2 photo-electrons 
pre  = precipitated  electrons 
gcr  = galactic  cosmic  rays 


The  ionization  due  to  energetic  particles  was  computed  by  first 
adding  the  pair  production  rates  of  the  source  models  of  galactic  cosmic 
rays,  photo  electrons,  and  precipitated  electrons  and  then  dividing  this 
pair  production  among  the  atmospheric  species.  Following  Swider  we 
assumed  that  the  cross  section  of  an  oxygen  atom  was  1.15  times  that 
of  a nitrogen  atom,  that  the  cross  sections  of  molecules  were  stoichio- 
metric multiples  of  the  cross  sections  of  their  constituent  atoms,  and 
that  these  ratios  were  approximately  valid  for  all  three  sources.  We 
also  chose  to  ignore  all  atomic  species  except  Nitrogen  and  Oxygen. 

Then 

and 


where 

gcr  = galactic  cosmic  rays 

phe  = photo-electrons  and  multiply  charged  ions 

prc  = precipitating  electrons 

R'  = ionization  rate  constant  for  species  n by 
particles 


Pairs,.  . , = Pairs  + Pairs  . + Pairs 

total  gcr  phe  pre 


r.  = <N>  + 1.15  <0> 
n [N]  +1.15  [0]  Total 


Pairstotal  = total  pair  production  rate 
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< > indicates  the  number  of  N and  0 atoms  in  species  n. 


[ ] indicates  the  total  local  atomic  number  density 
of  N or  0. 

1 9 

Dyadichov  and  Kozlov  give  values  for  relative  production  rates  and 
branching  ratios  which  corroborate  those  given  by  Swider.  The  major 
difference  is  in  the  relative  cross  sections  for  photoelectron  ionization. 
Uyadichov  and  Kozlov  give  cross  sectional  ratios  of  0:  0?  : N.,  = 0.56  : 1.5 
: 1.15,  while  Swider  gives  0 : 0.,  : ^ = 0.575  : 1.15  : T.0,  the  ratios 

used  in  these  studies. 

19 

Dyadichov  and  Kozlov  also  indicate  that  the  photoelectron  ioni- 
zation of  N.,  is  accompanied  by  dissociation  reactions  with  rates  amount- 
ing to  40%  - 90%  of  the  ionization  rate.  They  further  indicate  that  all 
of  the  N produced  by  dissociative  ionization  and  75%  of  the  N produced 
by  dissociation  by  photoelectrons  is  N(  D).  In  the  code,  the  computed 
va^ue  for  N ("D)  is  expected  to  be  low  because  these  reactions  forming 
N("D)  were  not  included. 

Dissociative  Ionization  wjis  included  only  for  the  two  species  N-, 
and  0-,,  whose  branching  ratios  are  0.19  and  0.29,  respectively.  Similar 
branching  ratios  were  found  by  Dyadichov  and  Kozlov. 

Predissociation  is  a reaction  path  which  includes  a molecule 
absorbing  energy,  being  raised  to  an  excited  state,  switching  to  another 
excited  state,  and  then  dissociating.  This  takes  time, and  radiation  or 
collisional  quenching  are  possible  alternate  paths.  The  branching  ratio 
between  radiation  and  predissociation  is  constant,  but  the  collisional 
quenching  rate  is  a function  of  M,  the  local  number  density,  and  hence  is 
a strong  function  of  altitude.  Ideally  this  would  be  modelled  by  includ- 
ing a metastable  species  and  the  several  decay  paths.  Because  the 
reaction  set  was  already  large,  the  rates  for  the  alternate  paths  were 
not  well  known,  and  these  paths  were  only  important  for  their  effect  on 
the  dissociation  path,  we  modelled  predissociation  by  chosing  a quantum 
efficiency,  <j>,  which  seemed  reasonable  over  the  altitude  $gnge  of  interest. 
For  N,,  we  chose  a if  of  zero  because  the  rate  is  very  slow^-^and  is  only 
effective  at  high  altitudes.  For  NO  we  chose  a <p  of  uni^“  for  our 
bands  < 28,  X < 190  nm.  For  0?  we  chose  a <j>  of  one  half  for  our  bands 
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V.  N.  Dyadichov  and  S.  I.  Kozlov,  "Minor  Ion  Components  in  the  Dis- 
turbed Ionosphere  in  the  Altitude  Range  100-200  km,"  Cosmic  Res.  15 , 

No.  2,  218-218  (1975). 
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P.  M.  Banks  and  G.  Kockarts,  Aeronomy , Part  A,  Academic  Press,  N.Y., 

1973,  p.  176  and  p.  318. 
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S.  Cieslik  and  M.  Nicolet,  "The  Aeronomic  Dissociation  of  Nitric 
Oxide,"  Planet.  Space  Sci.  21_,  925-938  (1973). 
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M.  Ackerman,  F.  Biaumd  and  G.  Kockarts,  "Absorption  Cross  Sections 
of  the  Schuman-Runge  Bands  of  Molecular  Oxygen,"  Planet.  Space  Sci. 

_18,  1639-1651  (1970). 
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27-29  (167-190  nm)  and  a of  zero  for  our  bands  > 29  (A  > 190  nm) 
which  qgrresponds  to  predissociation  for  v'  > 3 in  the  Schuman-Runge 
bands.  For  our  bands  < 27  the  value  of  <J>  was  absorbed  into  o.  The 
resulting  0,  dissociation  had  an  altit^e  variation  which  agreed  well 
with  that  given  by  Banks  and  Kockarts. 


PHOTOLYTIC  CROSS  SUCTIONS 


The  cross  sections  used  in  computing  the  photolytic  reaction  rate 
constants  were  chosen  at  times  when  our  knowledge  of  negative  ion  cross 
sections  was  changing  rapidly  and  hindsight  indicates  that  some  guesses 
and  extrapolations  from  private  communications  were  more  fortunate  than 
others.  Several  of  the  less  fortunate  choices  and  some  other  problems 
which  arose  in  choosing  cross  sections  will  now  b^  discussed.  Reaction 
numbers  without  parentheses  refer  to  the  ECLIPSE  reaction  set  and  will 
be  found  with  their  references  in  Tables  4 and  5 of  th^  Appendix.  Reac- 
tion numbers  within  parentheses  refer  to  the  BENCHMARK  reaction  set  and 
will  be  found  with  their  references  in  Tables  7 and  8 of  the  Appendix. 
The  actual  rate  constants  computed  for  BENCHMARK  are  shown  in  Table  9 
of  the  Appendix. 


For  the  reaction: 


NO- 


NO^  + e 


3 (1) 
-1 


a zero  shift  seems  to  have  occurred  in  the  conversion  from  cm  ' to 
band  number.  The  cross  sections  are  correct  for  the  next,  numerically 
higher  band  number.  In  effect,  the  cutoff  has  been  shifted  to  a 10% 
shorter  wavelength.  The  solar  flux  increases  strongly  in  the  cutoff 
region,  so  that  the  reaction  rate  constant  was  reduced  by  one  third,  and 
the  rate  constants  for  reaction  (1)3,  in  Table  9 should  be  increased  by 
50%  to  obtain  the  correct^^alue.  We  based  our  cross  sections  on  those 
computed  by  Herbst  et  al.  for  a vibrational  temperature  of  250  K.  We 
did  not  use  their  measured  cross  sections  because  thq|e  were  estimated 
to  correspond  to  an  effective  temperature  of  1750  K.“  Use  of  the 
measured  (1750  K)  cross  sections  would  increase  the  rate  constants  for 
Reaction  3(1)  by  an  additional  20%.  Herbst  et  al.  did  not  include 
rotation  effects  in  their  cross  section  calculations.  This  neglect  of 
rotation  should  be  equivalent  to  lowering  the  effective  temperature  by 
an  insignificant  amount. 
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Herbst  et  al.  also  measured  a photodissociation  cross  section  for 
the  peroxy  isomer  of  NC^  . This  species  is  not  in  our  reaction  set,  but 

23P.  M.  Banks  and  G.  Kockarts,  Aeronomy,  Part  A,  167-174,  Academic 
Press,  N.Y. , 1973. 
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E.  Herbst,  T*  A.  Patterson,  and  W.  C.  Lineberger,  "Laser  Photodetach- 
ment of  NO  “ ,"  J.  Chem.  Phys.  61_,  No.  5,  1300-1304  (1974). 

+ 1 

Read  as  Reaction  3 of  ECLIPSE  studies  and  Reaction  1 of  BENCHMARK  studies. 
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•>5  ">6 

could  be  important"  *"  because  the  photodissociation  rate  for 


N02  + hv  NO  + 0 


3a  (la) 


would  be  twice  the  corrected  rate  constant  for  Reaction  3(1).  One  possible 

1 * ‘ l_  C VIA  " *_  . 1 .• 


production  path  for  NO,  is  the  reaction 


C03  + NO  ->■  NO,  + C02. 


For  the  reaction 


co2  + O', 


later  measurements  and  estimates  by  Peterson  give  values  of  o which 
vary  from  135%  of  our  a at  longer  wavelengths  to  0.0%  at  shorter  wave- 
lengths. The  net  effect  is  that  our  computed  rate  constants  are  7-9% 
higher  than  they  would  have  been  if  we  had  used  Peterson* sZ  cross 


higher  th 
sections. 


For  the  reaction 


C03'(H,0)  + hv 


+ h2o. 


-17  2 

we  used  a constant  10  megabarn  (10  cm  ) cross2gection  throughout  the 
visible  ultraviolet  wavelength  range.  Peterson  gives  cross  sections 
which  are  lower  by  50%  in  the  longe|t  wavelength  band,  lower  by  75-80% 
between  400  and  520  nm,  and  zero  cm  for  X < 390  nm.  From  Peterson's 
cross  sections  we  would  compute  a rate  constant  for  reaction  8(7)  which 
is  about  45%  lower  than  the  rate  constants  which  we  obtained  from  the 
cross  sections  we  did  use. 

For  the  reaction 


N03  + hv  -*■  NO  + 02 

we  used  the  same  shape  at  cutoff  as  for  the  reaction 


12(12) 


NOj  + hv  -*  N02  + 0 + e.  4(11) 

Because  the  solar  photon  flux  decreases  muejj  more  rapidly  at  the  shorter 
wavelengths  where  the  N03  , 4(11),  and  N03~  , 12(12),  photodissociation 
reactions  cut  off  th^n  it  does  at  the  longer  wavelengths  where  the 
NO,  , 3(1),  and  N02  , 3a(la) , reactions  cut  off,  the  difference  between 

25 
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n 


the  rate  constants  is  much  greater  for  the  former  reactions  than  for  the 
latter  reactions. 

The  effects  of  these  changes  to  the  rate  constants  in  Table  9 of 
the  Appendix  arc  shown  in  Table  2. 

Positive  Ion  Reactions 


Positive  ion  photolytic  reactions  were  formally  included  in  the 
BENCHMARK-7b  reaction  set  in  Table  9 of  the  Appendix  as  reactions  18- 32 
but  in  practice  the  rate  constants  were  set  to  zero.  They  were  not 
included  in  the  ECLIPSE-66  reaction  set. 


Neutral  Species  Reactions 


The  square  root  law  variation  in  the  absorption  cross  section  of  0 ^ 
which  is  discussed  under  attenuation  computations  is  also  reflected  in 
the  07  dissociation  cross  sections.  For  normal,  Beer’s  law  absorption 


and 


I 


dl/dN 


I 

o 


no 

e 


- I 

o 


ae 


-no 


where  o is  a constant.  But,  for  square  root  law  absorption 

o = oq  (Nq/N) 1/2 
and 

dl/dN  = - 1/2  Ioae"n° 


so  that  the  effective  cross  section  for  the  dissociation  calculations  is 
1/2  of  the  effective  cross  section  for  the  attenuation  calculations.  Kith 
this  change  the  expression  for  the  effective  cross  section  for  07  disso- 
ciation in  the  square  root  law  region  becomes 


a 

e 


1.4 


a /(Jin (a  /a  . )(1 
max  v max  min 


2N 


a )1/2) 

max 


+ a . 
min 


for  all  N. 


Again,  due  to  an  algebraic  error,  the  factor  2.0  rather  than  1.4 
was  used  in  the  actual  computations.  For  0^  this  would  tend  to  compensate 
for  an  identical  error  in  the  attenuation  calculations.  The  absorption 
cross  section  error  raises  the  altitude  at  which  absorption  occurs  by  2-4 
kilometers.  The  dissociation  cross  section  error  raises  the  contribution 
of  the  affected  bands  by  50%.  For  09  the  latter  effect  would  predominate 
above  the  altitude  of  unit  optical  depth  and  the  former  effect  below  that 
altitude. 


25 


4 


Effects  similar  to  the  interaction  between  the  0,  absorption  cross 
section  and  the  CO,  absorption  cross  section  in  the  Sclmman-Rungc  region, 
our  solar  bands  27-31  (107-210  nm) , should  now  be  considered.  Both  CO 
and  11,0  have  photodissociation  cross  sections  which  decrease  very  strongly 
with  Increasing  wavelength  in  this  region.  Therefore  for  the  two  photo- 
dissociation  reactions 


CO,  + hv  ->  CO  + 0 16  (-) 

11,0  + hv  ->11  + 110  30  (34) 

we  made  adjustments  in  the  photodissociation  cross  section  identical  to 
the  adjustment  in  the  CO,  absorption  cross  section  which  was  discussed 
under  atmospheric  attenuation  calculations. 

The  solar  flux  is  increasing  moderately  with  increasing  wavelength 
at  the  same  time  that  the  0,  absorption  cross  section  is  decreasing,  an 
effect  which  partially  offsets  the  cross  section  interactions.  A correc- 
tion similar  to  the  cross  seetjon  correction  was  applied  to  the  photon 
flux  for  the  ECLIPSE  66  study.  It  overcompensated  at  low  absorption 
and  was  not  used  for  BENCHMARK  76.“ 

The  X-Ray-  Ionization  Cross  Sections  for  the  species  N,  N,,  NO,  NO,, 

0,  0,  and  0,(*A)  were  computed  from  the  N^  and  0,  X-ray  absorption  cross 
sections  using  the  formula 


where 


a.  . 
ij 


.5  <N>o..  . + .5  <0 >o„  . 

n2j  0,j 


a is  the  cross  section 


i is  the  reaction  number 
j is  the  solar  band  number 

< > refers  to  the  number  of  atoms  of  N or  0 in  the  species. 

The  9^os^  section  for  ionizatioq,of ,0,  by  solar  Lyman  B was  set  at 
1.6  x 10  “ m“  rather  than  1.0  x 10  m . This  would  have  a maximum 
effect  at  100  km  where  the  total  solar  ionization  rate  was  25%  high,  but 
the  effect  was  less  than  1/2%  at  90  km  and  negligible  below  that  height. 


This  reaction  was  not  in  the  BENCHMARK-76  set  because  CO  was  not  included 
in  the  species  set. 
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COMMENTS  ON  Till-:  CALCULATIONS 


Atmospheric  chemistry  models  are  frequently  allowed  to  run  for 
severaj ^mode 1 weeks  to  eliminate  all  start  up  transients.  The  A1RCI1LM 
models  were  designed  to  study  the  reaction  of  the  atmosphere  to  ioniza- 
tion shocks.  They  are  zero  dimensional,  transport  less,  and  make  no  claim 
to  model  slow  chemistry.  For  the  ECLIPSE-661  study  we  turned  the  sun  on 
at  10  AM  model  time  and  let  the  program  run  till  2 PM  model  time.  The 
2 PM  ion  densities  were  used  as  initial  densities  for  a second  run 
starting  at  10  AM.  The  new  2 PM  ion  densities  were  then  used  as  initial 
densities  for  a third  run  with  the  ECLIPSE  turned  on.  The  10:30  AM  model 
densities  usually  agreed  to  better  than  2-5%  and  the  2 PM  densities  for 
the  first  and  second  days  agreed  to  a few  tenths  of  a percent.  The 
neqtral  densities  stayed  within  their  estimated  error  bars  except  for 
N(~D).  This  density  fell  to  zero  because  several  formation  reactions^ 
were  not  considered.  (See  above  for  findings  of  Dyadichov  and  Kozlov.  ) 
For  the  ECLIPSE-66  study  the  SOURCE  package  was  called  to  update  the 
reaction  rate  cyijistants  at  both  the  prediction  and  correction  stages  in 
the  integrator.  For  the  BENCHMARK-^  study”  the  package  stood  alone, 
and  the  initial  neutral  density  model  was  used  to  compute  the  rate 
constants,  which  are  given  in  Table  9 of  the  Appendix,  and  which  also 
appears  as  Table  3 of  reference7^(I) . The  atmospheric  temperature  model 
was  taken  from  the  WEPH  4 code. 


-^W.  S.  Knapp,  R.  L.  Bogusch,  M.  I.  Chock,  "A  FORTRAN  Code  for  the  Cal- 
culation of  Ionization  and  Absorption  Due  to  a Nuclear  Detonation 
(WEPH  4),"  GE-TEMPO,  68  TMP-44,  DASA  2096,  30  April  1968. 
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APPENDIX 


This  appendix  contains  several  tables  of  input  data  used  in  the 
ECLIPSE-66  and  BENCHMARK-76  studies,  and  also,  as  a sample  output  the 
table  of  quiescent  production  rates  used  in  the  BENCHMARK  study.  Tables 
A1  through  A3  apply  to  both  studies;  the  other  tables  apply  to  only  one 
study  as  noted. 

Table  A1  lists  for  each  band  in  the  solar  model,  the  band  limits  in 
_g 

millimicrons  (10  m) , the  absorption  cross  sections  of  the  absorbing 

-2  - 1 

species,  the  solar  flux  in  photons  cm  sec  , and  the  model  band  number. 

Table  A2  lists  for  bands  11-20  the  ultra-violet  ionization  cross 

~22  2 

sections  in  megabarns  (10  M ) for  the  species  NO,  NC^,  09,  O.^A) 
versus  band  number. 

4 

Table  A3  lists  the  electron  precipitation  model  in  ion  pairs 
cm  3 s 1 versus  altitude  in  km. 

Table  A4  lists  for  bands  11-47  the  photodetachment  and  photodissocia- 

_2 

tion  cross  sections  in  cm  for  reactions  1-37  in  the  ECLIPSE-66  model. 

Table  A5  lists  for  all  bands  and  photolytic  reactions  the  references 
for  the  ECLIPSE-66  model. 

4 

Table  A6  lists  the  uv  and  x-ray  eclipse  functions  used  in  the  ECLIPSE- 
66  model. 

Table  A7  lists  for  bands  11-47  the  photodetachment  and  photodissocia- 
tion cross  sections  in  cm'^  for  reactions  1-49  in  the  BENCHMARK-76  model. 

Table  A8  lists  for  all  bands  and  photolytic  reactions  the  references 
for  the  BENCHMARK-76  model. 

Table  A9  lists  the  quiescent  production  rates  derived  from  the  above 
mentioned  data  for  the  BENCHMARK-76  model. 
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TABLi:  A2 . ULTRAV  IOLI'T  IONIZATION  CROSS  SECTIONS  EOR  ECLIPSU-66 
AND  BENCHMARK- 76 


SPECIES 

NO 

NO? 

°2 

o2(‘a) 

ECLIPSE  REACTION  NO. 

39 

40 

42 

43 

BENCHMARK  REACTION  NO. 

51 

52 

57 

60 

BAND  NO. 

11 

14. 

3. 

1.6* 

5. 

12 

8.5 

3. 

0. 

5. 

13 

6. 

3. 

0. 

4. 

14 

4.8 

3. 

0. 

3. 

IS 

3.2 

1.6 

0. 

0.7 

16 

2.4 

0.4 

0. 

0. 

17 

2.1 

0.12 

0. 

0. 

18 

1.9 

0.06 

0. 

0. 

19 

1.3 

0.004 

0. 

0. 

20 

0.7 

0. 

0. 

0. 

* 

This  was  apparently  a 

misprint 

which  occurred  when  the 

program  was 

punched.  1.0  x 10“  meter s^  is 

the  value 

indicated  by 

the  referenc 
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TABI.i;  A3.  ELECTRON  PRECIPITATION  MODEL  FOR  ECLIPSE-66  AND  BENCHMARK- 76 


Altitude  Pair  Production  Rate 

(km)  (ion-pairs  cm  ^ s *) 

100  .0047 
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TABLE  At>.  ULTRAVIOLET  AND  X-RAY  ECLIPSE  FUNCTIONS  FOR  ECLIPSE -66 


TIME 

UV 

X. 

X 

in 

out 

0 

.000 

.170 

. 165 

100 

.008 

.175 

.165 

200 

.018 

.180 

.170 

300 

.033 

.190 

.170 

400 

.044 

.200 

.175 

500 

.068 

.210 

.185 

600 

.089 

.225 

.190 

700 

.112 

.240 

.200 

800 

.136 

.255 

.215 

900 

.161 

.270 

.225 

1000 

.188 

.290 

.240 

1100 

.215 

.315 

.250 

1200 

.243 

.340 

.265 

1300 

.272 

.370 

.285 

1400 

.302 

.405 

.300 

1500 

.331 

.455 

.320 

1600 

.360 

.505 

.340 

1700 

.390 

.545 

.365 

1800 

.420 

.575 

.380 

1900 

.449 

.600 

.395 

2000 

.479 

.630 

.410 

2100 

.509 

.650 

.430 

2200 

.538 

.670 

.445 

2300 

.568 

.685 

.460 

2400 

.598 

.700 

.470 

2500 

.628 

.720 

.485 

2600 

.656 

.735 

.495 

2700 

.683 

.750 

.505 

2800 

.709 

.770 

.520 

2900 

.735 

.780 

.535 

3000 

.760 

.790 

.555 

3100 

.784 

.805 

.585 

3200 

.807 

.820 

.625 

3300 

.828 

.830 

.700 

3400 

.849 

.840 

.745 

3500 

.869 

.855 

.785 

3600 

.888 

.865 

.810 

3700 

.906 

.875 

.830 

3800 

.922 

.890 

.850 

3900 

.937 

.905 

.870 

4000 

.950 

.925 

.880 

4100 

.962 

.955 

.905 

4300 

.881 

.965 

.915 

4400 

.988 

.975 

.925 

4500 

.993 

.985 

.935 

4600 

.997 

.995 

.940 

4700 

.999 

1.0 

.945 

4800 

1.0 

1.0 

.950 

Time  is  time  to  or  time  from  totality  as  appropriate. 
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